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diferencia de potencial medida se obtiene por superposición de los efectos de las distintas 
fuentes o electrodos de corriente. 


Para la solución del problema fundamental se han aplicado básicamente dos procedi- 
mientos. El primero de ellos consiste en utilizar el método de imágenes, desarrollado por 
Maxwell [4]. Este método es de fácil aplicación al caso de un terreno formado por sólo 2 estra- 
tos. Sin embargo, se torna más complejo su uso para un número mayor de estratos. No 
obstante, se encuentran en la literatura ejemplos de aplicación para 3 y 4 capas [18, 211. 


El segundo método, utilizado 
con preferencia, consiste en la inte- 
gración de la ecuación de Laplace, 
sujeta a las condiciones de contorno 
presentes en un sistema estratifi- 
cado como el descrito. Este método 
fue aplicado por primera vez por 
S. Stefanesco y- los hermanos 
C. y M. Schlumberger en 1930 [16]. 


De acuerdo a este método, la 
solución general del potencial en un 
punto ubicado a una distancia radial 


r de una fuente puntual de corriente 
P, I, es: 


Figura 7.31 e Parámetros de un terreno estratificado. 


pi! l T 
A= elpt 20, RR | (7.29) 


donde: 
- 04 ¡(K, E, Y es la función kernel o función característica de Stefanesco, determinada por las 


condiciones de borde en las interfases entre los distintos n estratos que conforman el terreno. 
El primer subíndice se refiere al primer estrato. En lo sucesivo se suprime este subíndice, ya 
que para las mediciones de resistividad sólo interesa el potencial en la interfase aire-superficie 


del terreno. 
- Jo(r A) es la función de Bessel de primera especie y orden cero. 
- À es la variable de integración. 
- p4 es la resistividad del estrato superior del terreno. 
- I es la corriente inyectada en el origen O. 


La función kernel 6, (K, E, 4) es función de la variable de integración A, de los coeficientes de 
reflexión en las interfases entre estratos: 


Piel — Pi 
JS 7.30 
P a vaa 


y de los espesores E; de los estratos (1 < i < n); y se puede expresar como: 
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n = KEU A ON (7.31) 


H,(u)- P (u) ' 
donde: u=? ý 


A su vez, P,, y H, son polinomi i 
¡Y An polinomios en u, que pued i resiones recur- 
sivas siguientes: que pueden obtenerse aplicando las expre 


Ry 


I 


P (u) + Hu") K u" 


Hia: = H EBUS K, (7.32) 


con P4 = 0, H4 = 1. 
Por otro lado: h; = E, + E, +++ E), es la profundidad hasta el final del estrato i. 
De este modo, para una estructura formada por 2 estratos: 

HSHP KETE] 


Por tanto, para este caso, la función característica de Stefanesco es: 


0, = a P (7.33) 
o NS E A 
En forma similar, para una estructura formada por 3 estratos, se obtiene: 
-2E -2M E+E) 
Kie + Ke Gase 
Oy = l 2 (7.34) 


Tes K eE ar K, e 2ME¡+Ez) + Ki K eE, 


La expresión general (7.29), consta de una primera parte, correspondiente a la solución 
del problema en un medio homogéneo: término 1/r dentro del paréntesis cuadrado, más un 
efecto perturbador, por el hecho que no se trata de un medio homogéneo, representado por la 
integral de Hankel del lado derecho de (7.29). 


De acuerdo a la identidad siguiente [24]: 


oo 


[97 e Jy(rA) dk = (2? + pa , para m =0 
0 
= r(2 od ,param=1 , (7.35) 
el término 1/r de la expresión (7.29) puede escribirse como: 
1 o9 
Zs Ly (7.36) 
0 


Reemplazando (7.36) en (7.29), se obtiene otra forma de la expresión general: 
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pF i pl i 
(r) = e [1.+20,)/0(12)42 = er l| Np JolrA)dà, (7.37) 


es la denominada función kernel o función característica de Slichter, la 


donde N, =1+20,,, AN 
que incluye la solución fundamental del problema. Esta forma de expresar la solución directa, 


tiene ventaja en la obtención de la solución inversa; esto es, partiendo de valores medidos de 
resistividad aparente de un terreno, determinar la función caracteristica N,,(A) y, de ésta, obtener 
la función de resistividad p(z). 
A modo de ejemplo, la función característica de Slichter para una estructura de 2 capas 
es: 
K, e PE 1+ K eE 


x Sapoe ele = 77 (7.38) 
N, =1+20,=1 -Kei Ke E 


l 


Para una estructura de n capas, la función caracteristica de Slichter puede calcularse 
también mediante un algoritmo propuesto por Sunde [1]. 

El caso particular de un sistema de 2 capas, puede resolverse desarrollando en serie la 
expresión (7.33). Simplificando la notación a: 0, = 0 , Ķ4 = K , E4 = h4 = h: 


20h z = 
Ke = Keh y ketr A Y 2 (739) 


Ed ri 


Utilizando la identidad (7.35), la integral de Hankel contenida en la expresión (7.29), 


resulta: 
fos À)dà = [Y Kreh y (rà)dà = y e (7.40) 
olr ) = 0 = a 5 : 
0 0 n=1 n=lalr + (2 n h) 
Por tanto: 
IJl + K” 
Pı (7.41) 


a A SÓ 
q 27 |r o 


Una solución similar a (7.41) se obtiene aplicando el método de las imágenes a un siste- 
ma de 2 capas. En [4] se realiza un desarrollo detallado de aplicación del método de las imáge- 


nes para una estructura de 2 capas. 


7.7.2.2 Resistividad aparente 


Sobre la base de la expresión (7.29), que establece el potencial producido por una fuente 
puntual de corriente en un punto sobre la superficie del terreno, es posible determinar expresio- 
nes generales que permiten calcular los valores de resistividad aparente que deberían obtenerse 


al realizar mediciones en una estructura determinada. 


7 «AR 
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- Configuración de Wenner 


En términos generales, para la configuración de Wenner, la resistividad aparente teórica 
que deberla obtenerse de las mediciones es: 


274 
Pam = — (pa = 42) = Tor = A) — o(r=24)] (7.42) 


Por lo tanto, reemplazando (7.29) en (7.41), se obtiene: 


rd 


(Pa /Pi)m = 1+44f8,(K, EA) IA) Jo(241)]2 (7.43) 
0 


) | racterísticas de la estructura de un terreno: número de estratos, sus 
resistividades y espesores, la expresión general anterior permite determinar los valores de resis- 
A entre electrodos. 


Si se trata de una estructura formada por 2 estratos, procediendo en forma similar a lo 


anterior, pero utilizando (7.41) en lugar de (7.29), se obtiene la siguiente expresión para la 
resistividad aparente: 


ea K” 
(Pa! Pilian = 1+ 4) HH - 


K” 
1+(2nhl A)}  \4+(2nhl AY 


(7.44) 


- Configuración de Schlumberger 


Si se supone una separación pequeña entre los electrodos de potencial, la resistividad 
aparente teórica, según lo indicado en 7.3.4.2, es: 


L 
Pals) = PRE i (7.45) 
donde: E-=:=2(0010r). q, (7.46) 


es el campo eléctrico en el centro de medición. El factor 2 representa el efecto de ambas 
fuentes de corriente (electrodos C4 y C3). 


Introduciendo (7.29) y (7.45) en (7.44): 


2. d Jo(rh 
Pas = C, FE: ajo MAA a : (7.47) 
| 0 r=L 
donde: a = -AJ (rÀ) 
E 


Por tanto, para la configuración de Schlumberger, la resistividad aparente en relación con 
la del primer estrato, en función de la distancia L al centro de medición, es : 
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159) 


(Pa [Piho = 1420 20,0) 41090), (7.48) 


0 
donde J¡(LA) es la función de Bessel de primera especie y orden 1. 


Para el caso particular de una estructura formada por 2 estratos, sustituyendo (7.39) en 
(7.48) y aplicando (7.35): 


oo yr 


(Pa /Pi)s = 1+ 2% 


SEEEN EES (7.49) 
3/2 
n=1 ( + (2nh / Ly) 


En ciertas situaciones, al utilizar la configuración de Schlumberger, es posible que no se 
pueda cumplir con la condición de utilizar una separación a que sea pequeña en comparación 
con L. Para estos casos, los valores teóricos de resistividad aparente se determinan con (7.18), 
donde la diferencia de potencial Y =2[9(1 = L-a/2)- b(r= L+a/2)), se obtiene con 
(7.29) para la situación general, o con (7.41) para un sistema de 2 estratos. 


En [4] se analiza la diferencia entre los valores de resistividad aparente determinados con 
valores de a que no son suficientemente pequeños en comparación con L, y los valores límite 
teóricos dados por las expresiones (7.48) y (7.49) y con los cuales se construyen los gráficos 
patrón. Sobre esta base, es posible realizar una corrección aproximada de los valores medidos 
de resistividad aparente para la configuración de Schlumberger, para valores pequeños de L: 


4 
Pa(s)corregido = Pa(S)medido 4-(a/ L? 


7.7.2.3 Obtención de las curvas patrón de resistividad aparente 


Para la aplicación del método directo de interpretación de las medidas de resistividad 
aparente se requiere contar con las curvas patrón correspondientes, o disponer de un método 
que permita calcular los valores teóricos de resistividad aparente, si se aplica un procedimiento 
gráfico directamente en un computador. Si la estructura medida es de 2 capas, o se puede 
asimilar con aproximación aceptable para nuestros fines a esa condición; entonces los valores 
teóricos de resistividad aparente se pueden obtener de las expresiones (7.44) o (7.49), lo que no 
reviste mayor dificultad. Sin embargo, si se trata de un terreno formado por un número mayor de 
estratos, es necesario aplicar algún método que permita resolver las integrales de Hankel que 
aparecen en las expresiones de resistividad aparente antes vistas. 


En la literatura atinente se han propuesto diversas formas de solución del problema indi- 
cado [18, 21, 23]. Una manera consiste en aplicar el método de imágenes, lo que resulta labo- 
rioso para más de 2 estratos. Otra forma es efectuar una integración numérica de las expresio- 
nes (7.43) y (7.48), mediante Simpson o alguna variante, lo que se complica y produce errores, 
dado el carácter oscilatorio y poco amortiguado de las funciones de Bessel. Finalmente, otra 
opción es expandir la función característica (7.31), de manera que sea posible una integración 
directa de la serie resultante. 


Dentro del último grupo mencionado de posibles soluciones al problema, se describe 
brevemente a continuación el método propuesto por Mooney, Orellana y otros, consistente en 
expandir la función característica a un número suficientemente grande de términos, tal que se 
garantice una adecuada exactitud de los valores obtenidos [23, 5], j 
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De acuerdo al método mencionado, la función caracteristica de Stefanesco (7.31) 
expande como: 


9, = LON)“ (u= Ei (7.50) 


donde O(N) es un coeficiente a determinar, función de las resistividades de los estratos y SUS 
espesores, e independiente de la configuración utilizada y de las distancias 4, L. 


Sobre esta base, la expresión (7.29) para el potencial en un punto cualquiera sobre la 
superficie del terreno, se puede integrar directamente mediante las identidades (7.35): 


2 BI aa AN) 
0(r) = = + 2) === (7.51) 
2T? ME +(2NY 


De esta manera, las expresiones (7.43) y (7.48) de resistividad aparente toman la forma: 


(Pa Piw = = 1+4 (N) 


(7.52) 
N= Pra J+ 02N 1A? 
` 1 
(Pa Pi) = 1+ 2 Y O(N) (7.53) 


va (anti) 


Para calcular los coeficientes O(N), deben determinarse primero los valores de P(N), y 
H(N),, aplicando las relaciones recursivas siguientes: 


PON) = PON), FK HD- N); 
END) = AN); + K; P(D; - N); (7.54) 
con: P(0);= 0 , H(0);= 1, para todo i. 


Los valores de (7.54) se calculan hasta un valor de N =h „1, donde h „1 es la profundi- 
dad hasta el final del estrato n - 1. 


Una vez determinados los valores P(N)„ y H(N)y (0 <N < hp-1), los coeficientes Q(N) 
Se calculan con la expresión: 


d 
ON) = P(N) + LPU) =HON =a (7.55) 


j=l 
donde d es el menor de los dos enteros h„—] Y N-41. 
Lo anterior es sólo una descripción breve del método propuesto por Mooney, Orellana y 


Otros. Más antecedentes del método y de los criterios para su aplicación, pueden encontrarse en 
la bibliografía antes indicada. 
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7.4 métodos cuantitativos de interpretación 
7.7. 


Hasta la llegada masiva de los computadores, particularmente los de tipo personal, a ba 
la casi única alternativa de interpretación de las mediciones de resistividad aparente era la 
coso, ción de los gráficos de campo con curvas patrón. Sin embargo, en el área de la pros- 

e para usos geoflsicos, Se habían venido desarrollando, ya desde los años 30, numerosas 
e 


icas numéricas para una interpretación cuantitativa de las mediciones y obtener los paráme- 
En de la estructura investigada [22, 25 - 29]. 


Entre estas técnicas de interpretación, una de las más sencillas es aplicar el método de 
vadrados minimos al conjunto de mediciones de campo de resistividad aparente. Partiendo 
ls aa estimación inicial de los parámetros del terreno, el procedimiento consiste en ir modifí- 


do adecuadamente Éstos, hasta que los valores de resistividad aparente calculados se ajus- 
a io mejor posible a los valores medidos. 


Otros procedimientos más avanzados, aplican el método de los cuadrados minimos a la 
función característica de Slichter, descrita en 7.7.2.4, en lugar de los valores de resistividad 


aparente. La función característica se obtiene mediante una transformación numérica de las 
medidas de resistividad aparente. 


A continuación se muestra la aplicación del método de los cuadrados mínimos en el caso 
sencillo, pero práctico para nuestros propósitos, de una estructura formada por 2 estratos. 


tividad aparente. La estructura de terreno se supone constituida por dos estratos. Por tanto, se 


requiere determinar los parámetros Py, P2 Y E, = h, que correspondan a la mejor aproximación 
de las N mediciones. 


7 Pmed (A1), Pmed (Az) * + Ped (A y) son un conjunto de N mediciones de campo de resis- 


Para cada medición m Se define una función de error Fm: igual a la diferencia entre el 


logaritmo del valor calculado Pea y el logaritmo del valor medido Ped: Para UN determinado 


valor de la separación A, entre electrodos. La ponderación logarítmica de las magnitudes y 
parámetros proporciona una mejor solución del problema [22, 27, 29]. De este modo: 


F z LNP car Ar) = LNP mea (41) 
F, F LN Pear (42) Sx LNP med (42) 

z . i (7.57) 
cd LP An) LNP mea (An) 


Se desea determinar los parámetros Py, P2 Y h, que minimicen la función: 


N 


i | 
U) = Y [EF = YltnPealAn) — LP pea (An) (7.58) 


m=1 m=l 


Lnpy 
Y = Lnpa 


Lnh 


donde: (7.59) 


es la solución buscada. 


€ d 


Escaneado con CamScanner 


Pedro Ortuondo F MANUAL PARA EL PROYECTO Y ANALISIS DE PUESTAS A TIERRA 


Una forma posible de solución de (7.58) es utilizar el método de la mayor pendiente o 
método de la gradiente. Según [30], si Y” es una solución aproximada, una mejor estimación 


s A ad 
es. EA q A po 
ya! = yY? -pP WPR” (7.60) 
Fl 
p 
Lnp4” F2 
p =- 
donde: y? = |Lnp2” l F? =|. (7.61) 
Lnh? l 
p 
Fy 


(F? , WP WEF?) 


u’ = = (7.62) 
(W* yw. F’ y wW” w? ¿F?) 


donde, W?=W(Y?) y WP =W,(Y”) son respectivamente los valores numéricos del 


Jacobiano y del Jacobiano traspuesto, para y”. 


En las expresiones (7.60) y (7.62), el punto - denota producto matricial y (,) indica el 
producto escalar entre vectores. 


OF, 19Y, F/d% F 13V 
3F, 10%, ƏF,/Ə%, DF, 10% 


El Jacobiano es: W = ax. = (7.63) 
OFy 19Y, IFy10Y, ƏFy 10% 
donde, para cada fila m: 
DF, 19, = ILNP ea / 9, = Palð Peai / 9 P4) / Peat 
OF, / 9 Y == OLN Poy / DY, z P2 (0 Peat 1 0P2)/ Peal (7.64) 
ð F, /9Y, T IL Pe. / 0, sá HOP or 1 Oh) Poat 


Si en las mediciones de resistividad aparente se utiliza la configuración de Wenner, 
entonces, de acuerdo a (7,44): 


oo l 1 ` 
cal = l 4 > À "i —= - -= T 
Peal | + 2K E z) (7.65) 


donde: Q=1+ (nhl AY} (7.66) 
B = 4 + (2nh/ 4, y (7.67) 
K = (P3 =Pp1)/(p, + py) (7.68) 
7-60 
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y, por tanto: 


P) / h) ES e n l = K? 
Peal 101 = 1 +47, al si > es > $) (7.69) 
n=1 K a B 


dpal = 2 Yagi L 
Z lA (7.70) 


Pearl Oh = (16p4h I A7) Y $ K” n? tt (7.71) 
n=| i 


m 


Si en las mediciones de resitividad aparente se utiliza la configuración de Schlumberger, 
entonces, de acuerdo a (7.49): 


= l 
Pear = Py 1+2 K" -= ; (7.72) 
n=1 de 
donde: y =1+(2nh/L, y (7.73) 
y, por tanto: 
oo n 2 
dPear/ 3P = 1 + 2D, Sl ¿MENE (7.74) 
n=1 2K Ly 
E n(1-Ky? | 1 
Pear 1 9P2 = Y, a (7.75) 
n=l de 
agla = -Ap kI BIYAK" n ÁR (1.16) 


n=1 d5 


A partir de las expresiones anteriores es posible obtener iterativamente los parámetros p4, 
p2 y h de una estructura que se supone formada por 2 estratos. El proceso se puede detener 


una vez alcanzado un valor suficientemente pequeño de la suma U(Y), de los cuadrados de las 
diferencias entre los valores calculados y los valores leídos, según (7.58). 


La convergencia del método es buena y el número de iteraciones requerido para obtener 
una solución aceptable, depende lévemente de la estimación inicial de los parámetros del 
terreno, aunque esto no significa mayor problema considerando lo relativamente simple y rápido 
del proceso de cálculo. 
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Ejemplo 7.2 


Se desea determinar los parámetros de un terreno que se supone conformado por 2 estratos, 
cuyas valores de resistividad aparente, medidos con la configuración de Schlumberger, se 


indican en la tabla siguiente: 
m e p EEEE Dann, 
pogpgaga 


penie eeke 


Utilizando el procedimiento indicado anteriormente, en el gráfico siguiente se indican los valores 
medidos y la solución obtenida: 


p =2430m, Ej=h=43m, p,=380m 


| ZIF]? = 0.106D+00 


En rigor, los valores indicados de resistividad aparente corresponden con mejor aproximación a 
una estructura formada por 3 estratos, tipo K, con: p4 = 200 Q-m, E,=1m, p,= 400 Q-m, 
E, =2 m, y p3 = 40 2-m. Por tanto, los dos primeros estratos se han interpretado como uno de 
espesor E, = 4,3 m (mayor que E,+E>), y resistividad pı = 243 Q-m (menor que el valor medio 
ponderado de p, y p3). p> es con buena aproximación igual a la resistividad del tercer estrato. 


Una opción de interpretación de un sistema compuesto por más de 2 estratos, en particu- 
lar tipo H, K y similares, es utilizar un procedimiento por tanteo. Para esto se puede comparar en 
pantalla de un computador el conjunto de valores medidos, con curvas teóricas de resistividad 


+ A A 
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